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Abstract 

Anomalous interactions involving the top qnark contribute to some of the most 
difficult observables to directly access experimentally. They can give however 
a sizeable correction to very precisely measured observables at the loop level. 
Using a model-independent effective Lagrangian approach, we present the leading 
indirect constraints on dimension-six effective operators involving the top qnark 
from electroweak precision data. They represent the most stringent constraints 
on these interactions, some of which may be directly testable in future colliders. 


1 Introduction 

Once the presence of the Higgs boson has been firmly established at the Large Hadron 
Collider (LHC) the focus has turned towards the discovery of new physics (NP). Al¬ 
ready at Run 1 the absence of signihcant deviations from the Standard Model (SM) 
predictions has put stringent bounds on the mass scale of NP. Indirect constraints on 
new particles, when they are beyond the kinematic reach of the LHC, are becoming 
competitive and in some cases complementary to the indirect constraints from elec¬ 
troweak precision data (EWPD) pQ. It is expected that, with the increased energy 
available during Run 2, any new particle within the kinematic reach will be discovered, 
or very stringent constraints will be placed in case they cannot be directly produced. 
Still, there are certain interactions that, even with the increased energy, will be very 
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difficult to directly probe at the LHC with a significant precision. A notable example 
is that of the interactions involving the top quark that we consider in this work. 

Many well-motivated models of NP predict the largest deviations from the SM 
results in processes involving the top quark. Using an effective Lagrangian description, 
several groups have studied the potential of the LHC to constrain higher-dimensional 
operators involving the top quark in single and pair top production [2], including in 
some cases next-to-leading order predictions (see [3] and references therein). Some of 
these interactions can be directly constrained to a reasonable accuracy for the hrst time 
at the LHC. However, the complexity of the system limits the precision that one can 
achieve with these direct probes. In some cases, certain higher-dimensional operators 
are just inaccessible at the LHC. Many of these interactions, however, contribute at the 
loop level to EWPD. The very stringent constraints that can be derived from EWPD, 
together with the fact that the relevant coupling is usually the top Yukawa coupling, can 
compensate for the loop suppression, thus producing the most stringent constraints on 
many higher-dimensional operators involving the top quark. In particular, this includes 
those that cannot be directly probed at the LHC. 

In this article we use EWPD to place bounds on NP interactions. We use a 
model-independent effective Lagrangian approach going beyond the usual analysis of 
dimension-six operators correcting precision observables at the tree level. Analyses in¬ 
cluding one-loop contributions from higher-dimensional operators have been done for 
a small subset of purely bosonic operators in the past |1]. Here we use the calculation 
of the renormalization group equations (RGE) for the entire dimension-six effective 
Lagrangian [3 El E] (see also i) to determine, without restricting to any particular 
set of operators, which interactions can be constrained by EWPD at the 0(0.1) or 
better. This precision can be achieved by higher-dimensional operators that give a 
sizeable loop contribution to EWPD. If we further neglect operators that can be di¬ 
rectly probed in current or past experiments, we are basically left with dimension-six 
operators involving top quarks. We will show that indirect constraints from EWPD 
can be quite stringent for these interactions, some of which could be tested in future 
lepton colliders. 

We discuss the global fit to EWPD, including the leading, currently unconstrained, 
loop effects in section [21 The corresponding constraints on the coefficients of the 
dimension-six operators involving the top quark are presented in section [3] and we 
conclude in section 01 

2 The global electroweak fit for new physics to di¬ 
mension six 

Assuming that NP is heavier than the energies currently probed by experiments, its 
effects can be described by an effective Lagrangian, 

+ ^-^6 + • • • , ( 1 ) 


2 


where Cd = ^ CiOi, with Oi Lorentz and SM gauge invariant operators of mass dimen¬ 
sion d built from the SM helds. A stands for the cut-off scale of the effective theory, and 
we have neglected lepton number violating effects. There is a total of 59 operators (up 
to flavor indices) at dimension six [9]. However, only a few of these directly contribute 
at leading order to EWPD. We consider in our analysis the following set of EWPD; 


Mh, mt, as{Ml), Aa® (M|), Mw, Pw, Mz, Tz, cxhad, R 


Ae ,^,.( SLD ), 

Rb,c^ Ru/Ru-\-d-\-s^ Qf^| 


’e,/i,T5 

^had I 


( 2 ) 


whose dehnition, experimental values, errors and correlation matrices are taken from 

pro]. 

With these observables we perform a £t to the SM extended with the relevant set 
of dimension-six interactions that we will introduce below. The SM predictions for 
EWPD are computed including the latest theoretical developments HU- These are 
then corrected by the NP dimension-six operators, and compared to the experimental 
measurements via the usual function, which we minimize and use to determine 
the bounds on the new interactions. In our hts both the NP and the SM parameters 
are allowed to float. In order to include dimension-six effects in a consistent way, 
we limit the contributions to EWPD to order 1/A^. These may come either from 
the direct interference of the NP amplitudes with the SM ones in the electroweak 
precision observables, or from NP contributions to the physical processes where some 
of the SM input parameters are determined (e.g. Gf as extracted from muon decay), 
which then propagate to all observables. Effects of order 1/A^ would be comparable 
to interference effects coming from dimension-eight operators and have been therefore 
neglected. Given the resulting bounds, this seems in general a good approximation. An 
analysis of the error induced by neglecting such contributions, along the lines explained 
in e.g. ns, goes beyond the scope of this paper. 

The analysis of bounds on dimension-six interactions from EWPD has been con¬ 
sidered extensively in the literature (see for instance [I3]). At the tree level, only the 
following coefficients of dimension-six operators modify the observables in (l2i) : 


tree _ 
'-EWPD — 


db .(1) .(1) .(1) .(1) .(3) d3) 


( 1 ) 


^(pl 1 ^4>q I ^(pe 1 ^pui ^pdi ^pl 1 ^(pq’’ ^<pD) CWBj 


(3) 


The definition of the operator basis used in this work follows closely the original one in 
[9] except for the four-fermion sector for which we use the one in Appendix A of [TT] 0 
At the loop level, however, many other dimension-six operators contribute to EWPD, 
including some to which we have little or no direct experimental access. In this latter 
case, the high precision of the EWPD can compensate the loop suppression and give 
the most stringent constraints on these operators. The complete one-loop calculation 
of the EWPD including dimension-six operators is beyond the scope of this article 

^For completeness, the dimension-six operators in m are presented here: D^,(=((/)bD^b)(d’7^b); 
D0^,^=(dhD“(/))(bL7'^craV'L) (withZ?^ = D^-D^ andD^ = (JaD^-D^ao), O^D=i4>'^cj)), 
and (TTmG) 
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but the logarithmically enhanced contributions can be computed by means of the RGE 
recently computed in ISIEIIZ]- An analysis of the loop-improved electroweak constraints 
on the dimension-six interactions in ([3]) will be presented elsewhere [13]. Here we want 
to focus on those operators that, as explained above, have not been directly probed by 
any experiment yet, and to determine which ones can be constrained to a reasonable 
accuracy with current data. 


The largest RGE effects on EWPD are those proportional to the top Yukawa cou¬ 
pling, Ut. (The strong coupling does not enter in any of the anomalous dimensions 
for the interactions in dSj) .) Inspecting the RGE for the operators in ([3|) and looking for 
Ut effects then allows us to identify which ones could be significantly constrained from 
the low-energy bounds. Restricting to operators that have not been directly tested in 
experiments so far leaves us with the following set, containing always the top quarkH 


C’ig = i(p^'iD,.(p)iqLYqL), 

Ofh 

Oll^ = (k^,lL){qirqL), 

0g)= 

Oeu = {'^'li,eR){u]iYuR), 

Oqe =(^ 7 /. 9 L)(^ 7 ^eR), 

0,„ = 

= ^{qL'l^lqL){W'y^qL), 


C’iii = ^(Ml^,UR)(u]i-f>"uR), 

oS = 


07 = 

OuB = 

OuW = 




ilL7^l(TJL)iqL'J^(TaqL), 

^ (Ml^lTAqL) iWl^TAqi), 

i'^'y/^UR){dR'y^dR), 

{WltiqL){dR'y^dR), 

{Wcr^’'(^aUR)^W^^. 


(4) 


An analysis of EWPD constraints on a small subset of these operators was per¬ 
formed in HU- However, at one loop only non-logarithmic hnite contributions were 
included and the corresponding bounds are much weaker. Regarding the operators 
we will only consider the combination which (up to cor¬ 

rections suppressed by products of Vtd and Vts with V the Cabibbo-Kobayashi-Maskawa 
matrix) modifies the neutral current top couplings, without inducing any tree-level cor¬ 
rection in the bottom ones. As we will also see, because EWPD is only sensitive to 
= (C>^q^ -I- at the tree level, NP effects coming from Oqf cannot be con¬ 
strained by the data. Similarly, EWPD is not sensitive to so Ouu cannot be 

constrained if it only involves the third generation. 

As we mentioned above, we compute the predictions for physical observables consis¬ 
tently with the approximation of a dimension-six effective Lagrangian, including only 

^Apart from these, we will also comment on another operator (see footnote 3): the pure scalar 
interaction 4>), which does not contribute to EWPD proportionally to yt, but enters 

in the anomalous dimension of c^d with a large coefficient ~ 20gi/3 |7]. Although this operator can 
be in principle tested in Higgs physics, its constraints are still very weak m- 
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the interference between the SM amplitudes and the NP effects (i.e. terms linear in 
1/A^). We also use the leading logarithmic approximation for the solution of the RGE 


dCi 1 

d log fi 


tiCi Ci M « pf + 


Y^7/(A)log£)Q(A), (5) 


where we have dehned the dimension-six coefficients Ci = Ci/h?'. We include in the 
anomalous dimensions 7 ^ the full dependence on the SM gauge couplings [7j and 
the leading contributions from the Yukawa interactions [ 6 ], i.e. Yg,Yd ^ 0, and 
Yu ~ diag (0, 0, i/i). The dependence on the Higgs self-coupling [5] is irrelevant for 
our analysis. Note that, for the contributions to the EWPD that only arise from the 
RGE, within the approximation in Eq. ([5]) the physical predictions depend always on 
Cilogfi/A. Finite one-loop contributions beyond the logarithmically enhanced terms 
included in the RGE can in some cases be relevant [18], and therefore modify the pre¬ 
vious dependence. Whenever these hnite terms are available, we include them in our 
analysis. We explicitly comment on these cases below. 


3 Loop constraints on new top interactions 

In this section we present the constraints that EWPD impose on the operators involving 
the top quark in (jl]) due to RGE effects. We report these bounds in Table [U assuming 
that only one operator is generated at the ultraviolet scale. A, at a timejl We present 
the results from two different types of hts. The limits on CilogMz/A (and q, for 
A = 1 TeV) are obtained assuming that the dimension-six coefficients can have either 
sign. On the other hand, the bounds on the NP scale A are derived from a £t with 
the extra assumption that the c* have a dehnite sign, and then setting this to some 
illustrative values, Cj = ±1. As we mentioned above, the precision of EWPD overcomes 
the loop suppression of the RGE effects and allows to constrain most of the interactions 
at the few percent level for A = 1 TeV, or alternatively pushes the scale of NP in the top 
sector to a few TeV for q = ±1, thus fully justifying the use of an effective Lagrangian 
description. (Note that, even for the weakest bounds, the NP scale is always pushed 
signihcantly above the Z mass, where the EWPD are measured.) The leading-log 
approximation used is also justihed provided the value of A is not too large, so that 
I ^ log'll -C 1, with Q!j the relevant SM parameter. In the rest of this section we 
discuss the origin of the constraints on the different operators. 

The relatively strong constraints on interactions can be understood from 

the fact that all those interactions contribute to the running of or The 

corresponding operators provide direct corrections to the neutral current couplings of 

^For completeness, we include here the 95% probability interval for the operator coefficient c^n: 

^ log ^ e [-4.63,0.65] TeV^ (c^n S [-0.27,1.93] for A = 1 TeV). (6) 

For negative values of the coefficients, the corresponding bound is then somewhat better than the one 
obtained from Higgs observables [IT] . 
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Operator 

95% prob. 

interval 

95% prob. 

lower bound 



Ci 

A 

[TeV] 


[TeVf 

(A = 1 TeV) 

(Ci = +1) 

(Ci = -1) 


[-0.15,0.38] 

[-0.16,0.06] 

4.4 

3.2 

(olfUu 

[-0.26,0.36] 

[-0.15,0.11] 

3.7 

3.3 


[-0.21,0.44] 

[-0.18,0.09] 

3.8 

2.9 


[-0.40,0.16] 

[-0.07,0.17] 

3.1 

4.3 

(^ge)ttee 

[-0.42,0.20] 

[-0.08,0.18] 

3 

3.9 

(Of Am. 

[-0.91,0.25] 

[-0.11,0.38] 

1.9 

2.9 

(Of Am. 

[-0.04,0.54] 

[-0.22,0.02] 

4.8 

2.6 

eu) 

[-1.29,0.22] 

[-0.09,0.54] 

1.5 

2.6 


[-0.26,0.95] 

[-0.40,0.11] 

2.8 

1.9 


[-0.22,1.24] 

[-0.52,0.09] 

2.7 

1.6 

(Of Ax.. 

[-0.52,0.96] 

[-0.40,0.22] 

2.3 

1.8 

(Of Ax,. 

[-0.86,0.69] 

[-0.29,0.36] 

1.9 

2.1 

i^^eu^Trtt 

[-0.58,1.18] 

[-0.49, 0.24] 

2.1 

1.6 


[-1.01,0.54] 

[-0.23,0.42] 

1.8 

2.2 

(^^qejttrr 

[-1.14,0.56] 

[-0.23,0.48] 

1.7 

2.2 


[-0.16,0.26] 

[-0.11,0.07] 

4.7 

3.9 


[-0.07,0.29] 

[-0.12,0.03] 

5.9 

3.8 

( ) £itt 

[-0.24,0.33] 

[-0.14,0.10] 

3.8 

3.4 

{_^lu) ^£tt 

[-0.27,0.17] 

[-0.07,0.11] 

3.8 

4.6 

{Oqe)tm 

[-0.32,0.23] 

[-0.10,0.13] 

3.4 

3.9 

(of A... 

[-0.55,1.38] 

[-0.58,0.23] 

2.1 

1.5 

(ofA.» 

[0.25,10.9] 

[-4.6,-0.10] 

0.89 

0.37 

(ofA... 

[-1.47,0.59] 

[-0.25,0.62] 

1.4 

2 

(ofA,» 

[-9.7,-0.07] 

[0.03,4.06] 

0.41 

0.95 

{OuB)tt 

[-0.35,0.10] 

[-0.04,0.15] 

3.4 

5.1 


[-0.39,0.11] 

[-0.05,0.17] 

3.2 

4.7 


Table 1: EWPD bounds on top interactions, assuming one operator at a time at the 
scale A. The bounds on the NP scale A are obtained from two independent types of fits, 
assuming a dehnite sign for the coefficients q. The results for the operators 

are obtained assuming lepton universality in the interactions. The bounds for {Oqq )utt 
are too weak and have been omitted, while the operator coefficient for {Ouu)tttt cannot 
be constrained within our approximations (see text for details). 
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the charged leptons, and are bounded at the per mile level (see [19] for an earlier partial 
analysis). Note that the bounds for some of these interactions, e.g. and Oiu, are 
almost identical (up to a sign). This correlation follows directly from the RGE for the 
leptonic interactions in ([3]). In particular. 


d{c!;x 

dlogfi 

dlogfi 



(7) 


Thus, only the combinations of operators appearing in Eq. ([7j) can be constrained by 
EWPD, up to corrections in the RGE induced by the gauge interactions. 

Gonstraints on four-quark interactions involving only the third family are dominated 
by the contributions they generate to the Zbb couplings, via the operators and 

and are therefore somewhat weaker than the leptonic ones. Limits on c^q and 
Cqu arise from the bounds on the left-handed bottom couplings, and are signihcantly 
stronger than those of and , which contribute to the Zb^bji interactions. In 

particular, the strong preference for a positive (negative) value of {{c^qd)ttbb) 

follows from the corresponding preference for a large correction to the right-handed 
bottom coupling, = —^{c^^)bb^, with v ~ 246 GeV, to alleviate the —2.5-cr devi¬ 
ation in the bottom forward-backward asymmetry at the Z-pole. Again, some of the 
bounds on these four-quark operators can be easily correlated from the corresponding 
contributions to the running for the quark interactions in (jS]), 






+ ... , 


+ ..., 


( 8 ) 


that determine which combinations of operators can be constrained by EWPD. In 
the hrst line of Eq. ([8]) there is no contribution from Og^q , because the corresponding 
corrections to the running of and G® cancel each other. There is a suppressed 

contribution to the running of from Cqg^ proportional to the electroweak 

gauge couplings, which results in much weaker constraints. This explains the absence 
of a bound on Cqf in Table [H Finally, the coefficient {Ciru)tttt only contributes to the 
ZtRtR couplings through the RGE for and therefore cannot be bounded by 

EWPD at the order we are working. 

Four-quark operators involving two quarks of the third generation and two of either 
the hrst or second generations contribute, through RGE, to operators that modify the 
electroweak couplings of the quarks in the hrst two generations. These have been mea¬ 
sured with worse precision than those of the charged leptons or bottom quark. Hence, 
the corresponding bounds are much weaker and not reported here. If one assumes 
universality among the three families then the bounds are still mostly dominated by 


7 










the operators involving only third generation quarks. The exception is the case of 
the operators for which there is a tension between the required contribution to 

5g\i and the corresponding one for the hrst two generations. This tension results in 
signihcantly improved bounds, reducing the size of the corresponding 95% probability 
intervals by a factor of two, e.g. {C^qd )ttqq^og^ G [—4.09,0.65] TeV“^. 

The limits on the electroweak top dipole interactions, {OuB)tt and {Ouw)tt, come ex¬ 
clusively from their contributions to the running of ciys (cvi/s/A^ G [—0.009,0.003] TeV“^ 
at 95% probability), which is related to the S parameter [SU]. Hence, only the approx¬ 
imate combination g 2 {CuB)tt + 2(y'i(| + \){Cuw)tt (where the 1/6 and 2/3 factors are 
the qi and ur hypercharges, respectively), which enters in the RGE for cwbi can be 
constrained. 

Finally, we have not included in Table [T]the constraints on the operators that induce 
direct corrections to the top electroweak couplings. 


^91 = 


1 

2 






Hc'b — 
2 V <’“/« K‘' 


(9) 


Note that, to dimension six, the effects on the left-handed sector are also correlated with 
the direct corrections of the charged current couplings, = (Vc (We work 

in a flavor basis in which the SM Yukawa couplings for the down sector are diagonal.) 
The constraints on the combinations in Eq. ([9]) follow from the one-loop contributions 
to the T parameter and corrections to the Zbb vertices. These corrections contain 
logarithmically enhanced terms that can be read from the RGE of the operators 
(equivalent to the T parameter in our basis), + 0^^^)bb and {Ol^]^)bb- We 

have also included hnite (not proportional to logarithms) one-loop effects by integrating 
out the top quark with the anomalous couplings dehned in Eq. (P [2I1E2]. In particular, 
the hnite contribution to the T parameter is given by 



Because of these hnite terms, the depends on both G, and GRog^, so we vary 
Gj and A independently in our hts. We impose A > 1 TeV, to avoid regions where 
the ehective Lagrangian description may break down. In the bounds below, when no 
mention to A is made, we take the most conservative bound that is reached for A = 1 
TeV. 

Gonsidering only one of the combinations in Eq. ([1]) at a time we obtain the following 
95% probability interval for bg\/g\^^, 


AWI 0.016,0,002] 

9l 

(4‘>e [-0.01,0.09]), 

(11) 

while for bg^^/g^^^ we get 



^ 6 [-0.017,0.002] 

9r 

((A>)«e [-0.08,0.01]), 

(12) 
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t SM 




Figure 1: Left) 95% probability regions in the 

9l 9ji 

TeV. Right) Boundaries of the 95% probability intervals for 


plane for A =1, 3 and 10 

5g^ 

-Tw as a function of A. 

9l,r 


where ~ I “ | and = — | sin^ 9w, with 9w the weak angle. When we 

consider dimension-six effects correcting both electroweak top couplings at the same 
time, the 95% probability bounds change to: 

Sq^ S 

e 1-0.048,0.0891, ^ 6 1-0.102,0.044]. 

Sl Sr (13) 

(c“ e [-0.52,0.28], (c«)« e 1-0.50,0.21].) 

The weaker bounds compared with Eqs. fllll) and 0121 ) follow from a strong correlation 
(~ —99%) between these two couplings, as can be seen in Figure [H left, which can be 
understood from the leading logarithmic contributions. Indeed, neglecting the finite 
contributions, there is a large correlation between the corresponding dimension-six 
operators, which comes from the RGE for C^e> (whose limit dominates the constraints 
in Eqs. flTTD and flT^ via the T parameter). 


dC^D 

dlog/i 



(14) 


At this level there is an approximate flat direction for Sgf/g\f^ = —dg)i/gji^ (we have 
used gjf^ ~ —2gjf^ and = —c^^), which is also reflected in the equality of the 
bounds in Eqs. (na and flT^ . This flat direction is lifted however by the logarithmic 
contributions to the Zbibi vertex and, to a less extent, the contributions in the running 
from gauge interactions. The global factors in front of the finite terms turn out to be 
smaller than the ones in the logarithmic terms, so the effects from the former are not 
very important for values of A consistent with the effective Lagrangian description. In 
Figured] right, we show how the bounds obtained in Eq. flT^ for A = 1 TeV evolve as 
we increase the value of the scale of NP. 
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The results in Eq. imply quite strong bounds on deviations with the respect to 
the SM electroweak theory, at the few percent level. One may wonder though about the 
robustness of these bounds when interpreted within particular models. For instance, as 
explained above, the effects in the T parameter have a strong impact in these bounds. 
Now, the T parameter is known to have a strong correlation with the S parameter (we 
obtain an 89% correlation), and many SM extensions in which the NP interacts mainly 
with the third generation come also usually accompanied by relatively large positive 
contributions to S. Still, considering a positive value for S = 0.2 (about the 95% 
probability limit obtained from the S-T £t) does not have a dramatic impact on the 
bounds in Eq. flT^ . due to the constraints on the logarithmic contributions to Zbibi- 
We thus obtain only a moderate shift in the bounds for the right-handed couplings, 

^€[-0.050,0.0881, ^6 [-0.123,0.023] (S = 0.2), (15) 

i/L 

but we can still conclude that NP contributions to Zti couplings beyond 10% are 
disfavored by the data. 


4 Conclusions 

Among the possible NP deformations of the SM there are several on which no direct 
experimental information can be extracted from present or past experiments. Indirect 
constraints can however be obtained in some cases from their loop contribution to 
the very precisely measured EWPD. Using a model-independent effective Lagrangian 
approach, we have explored in this article the potential that EWPD have to put con¬ 
straints on dimension-six effective operators on which no precise direct information 
is currently available. We have shown that, despite the loop suppression, the lead¬ 
ing contributions proportional to the top Yukawa coupling are large enough to place 
signihcant bounds on several interactions involving the top quark. 

Our results, reported in Tabled! show that, for NP in the TeV region, EWPD can 
constrain the dimension-six operator coefficients for a large set of interactions involving 
the top quark to 0(0.1) values. These results are obtained using the logarithmically 
enhanced one-loop contributions, as given by the RGE, together with the hnite terms 
in the most signihcant cases. Barring accidental cancellations, similar bounds will 
apply for all the interactions we have considered, even if the missing finite terms give 
a comparable effect to the logarithmically enhanced ones. 

The bounds presented here have been computed assuming only one operator at 
a time at the scale of NP. There are several scenarios with new scalars, quarks or 
vector bosons for which such operators can be generated alone, or whose effects are 
not correlated with the contributions to other dimension-six interactions [Il|23l[2l]. 
Nevertheless, we have discussed the origin of the leading constraints and described 
the most relevant approximate hat directions so that bounds on more complicated 
models can be estimated. As an example, let us discuss the case of lepton-top four- 
fermion interactions. For NP at the TeV scale we are able to place 0(0.1) bounds on 
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the coefficients of contact interactions resulting from the products of vector currents 
involving two electron and two top helds, with all possible chiralities. As we have 
mentioned, the effects of some of these operators are correlated in the RGE that enter 
in the EWPD so, in practice, we can only constrain three combinations of these e^e~tt 
operators. Removing the two redundant interactions from a simultaneous global £t to 
all the operators we obtain the following limits (for A = 1 TeV): 


^Iq ^lu)eett 

e 

[-0.21,0.08], 


\^lq jeett 

e 

[-0.31,0.06], 


^eu)eett 

e 

[-0.10,0.29]. 

(16) 


In particular, these e+e tt contact interactions can be tested in future colliders 1251, 
like the ILC [26] or the FCC-ee [27| , and are therefore relevant to guide NP searches in 
these facilities. The same applies for our analysis of the Ztt couplings, which could be 
directly measured at this kind of experiments with great precision (around 1% for the 
ILC with 500 fb“^ |2B])- Finally, while less precise, our limits on four-top interactions 
are comparable (and in some cases signihcantly better) than the latest ones obtained 
by the LHC [29] . 

Higher-dimensional operators involving the top quark are not only among the least 
constrained operators that contribute the most to EWPD, they are also among the 
best motivated ones in models that attempt to solve the hierarchy problem. In many 
of these models, a sizeable contribution to the S parameter is generated at the NP scale. 
We have considered that possibility, see Eq. flT^ . and shown that while inducing a shift 
in the allowed values of the coefficients of the higher-dimensional operators involving 
the top, these can still be constrained to a similar level of accuracy. This also illus¬ 
trates the fact that, being indirect bounds, the limits we have computed are sensitive 
to assumptions about operators that can give tree-level contributions to EWPD. In the 
presence of such operators the quantitative results may change but the qualitative fact 
that one-loop contributions can place the most stringent bounds on currently untested 
operators still holds. These bounds are therefore a crucial piece of information, both 
for model building purposes and as a guide for future experimental searches. 


Note added: While this manuscript was being prepared for submission, ref. [30] 
appeared in the arXiv. In that work, a global fit to dimension-six effective operators 
involving the top quark is performed using top production data from the LHC and the 
Tevatron. Our approach is complementary to theirs, in the sense that our analysis is 
sensitive to a different class of operators. In the few cases in which there is overlap, 
our results provide more stringent bounds on the coefficients of the operators. 
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